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The Hohenberg-Kohn (HK) theorem is one of the most fundamental theorems of quantum mechan- 
ics, and constitutes the basis for the very successful density- functional approach to inhomogeneous 
interacting many-particle systems. Here we show that in formulations of density-functional theory 
(DFT) that employ more than one density variable, applied to systems with a degenerate ground 
state, there is a subtle loophole in the HK theorem, as all mappings between densities, wave functions 
and potentials can break down. Two weaker theorems which we prove here, the joint-degeneracy 
theorem and the internal- energy theorem, restore the internal, total and exchange-correlation en- 
ergy functionals to the extent needed in applications of DFT to atomic, molecular and solid-state 
physics and quantum chemistry. The joint-degeneracy theorem constrains the nature of possible 
degeneracies in general many-body systems. 



Introduction. Quantum mechanics is based on the as- 
sumption that all information that one can, in principle, 
extract from a system in a pure state at zero temperature 
is contained in its wave function. In nonrelativistic quan- 
tum mechanics the wave function obeys Schrodinger's 
equation, 1 which implies a powerful variational principle 
according to which the ground-state wave function min- 
imizes the expectation value of the Hamiltonian. This 
variational principle was used by Hohenberg and Kohn 
(HK) 2 to show that the entire information contained 
in the wave function is also contained in the system's 
ground-state particle density n(r). 

HK established the existence of two mappings, 



v(r) <^ #(ri,...rjv) 



n(r), 



(1) 



where the first guarantees that the single-particle poten- 
tial is a unique functional of the wave function, v[*&], and 
the second implies that the ground-state wave function 
is a unique functional of the ground-state density, 4"[n]. 
Taken together, both mappings are encapsulated in the 
single statement that the single-particle potential is a 
unique density functional v[n]. In this formulation, the 
HK theorem forms the basis of the spectacularly success- 
ful approach to many-body physics, electronic-structure 
theory and quantum chemistry that became known as 
density- functional theory (DFT)ii^ 

Mapping 2 was originally proven by contradiction^ and 
later by constrained search^ Note that, in spite of occa- 
sional statements to the contrary in the literature, nei- 
ther proof directly proves the combined mapping, and 
thus the existence of the functional v[n]. This requires 
additionally mapping 1 , which in the case of density-only 
DFT is proven by inverting Schrodinger's equatior^. 



According to the preceding equation, any eigenstate 
determines the single-particle potential v(r) up to an 



additive constant, the corresponding eigenenergy^ The 
ground-state version of this equation is all one needs to 
prove mapping 1, and thus of the combined mapping 
n(r) v(r) + const, which we write as v = v[n]. The 
constant can always be absorbed in the definition of the 
zero of energy, but the resulting nonuniqueness of the 
mapping v [n] gives rise to derivative discontinuities^ that 
crucially contribute to observables such as the band gap, 
chemical hardness, electron affinities, electron transfer 
energies, among others. 

The two classic proofs of the HK theorem, by contra- 
diction and by constrained search, have been extended 
to a wide variety of systems, including spin magnetiza- 
tion (in spin-density-functional theory - SDFTiSiAi) and 
orbital currents (in current-density-functional theory - 
CDFT 12 ), among others. Further scrutiny, however, has 
led to the discovery of situations in which each of the two 
mappings breaks down. 

Breakdown of mappings: Degeneracy. Mapping 2, 
*&[n], breaks down in the presence of degenerate ground 
states, where neither the constrained search proof nor 
the proof by contradiction provide a unique wave func- 
tion among the degenerate manifold, because the strict 
inequality of the underlying variational principle is re- 
placed by the weaker relation less-or-equal. In the case 
of the proof by contradiction, the contradiction simply 
does not follow unless one has a strict inequality, and in 
the proof by constrained search there is no guarantee that 
the search delivers only one wave function for a specified 
density. 

What remains from the proofs is that the ground-state 
density uniquely determines a manifold of degenerate 
states, {^^[n], but not all of these states individually. 13 
Any member of this manifold, however, still uniquely de- 
termines the potential, since any of them can be used in 
Eq. (0). Hence, even in the presence of degeneracy, the 
mapping v^i], and thus v[n], still exists. This situation 
is illustrated in black in Fig. ^ 

To proceed from establishing mappings to a practi- 
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cal density-functional theory, one must define the total- 
energy functional S„[n], the universal internal-energy 
functional F[n], and the exchange-correlation energy 
functional E xc [n] . The Kohn-Sham formulation addition- 
ally requires the nonintcracting kinetic-energy functional 
T s [n]. 21 years after proving the original HK theorem, 
Kohr. 14 showed how these functionals can be defined even 
if degeneracy renders the original proof ineffective. 4 Since 
all degenerate wave functions by definition yield the same 
ground-state energy E, one can directly define the func- 
tional 

F[n] := E - J d 3 rn(r)v[n]{r) = E - V[n). (3) 

Conventionally, this functional is defined as F[n] = 
T[n] + U[n], but the information that the kinetic energy 
T and interaction energy U are density functionals is only 
available if the second mapping, *&[n], holds, and cannot 
be taken for granted in the presence of degeneracy. By 
contrast, the alternative definition above only requires 
the mapping v[n] to establish the existence of the uni- 
versal internal-energy functional F[n\. Similarly, in the 
nonintcracting KS system one defines 

T s [n] := E s - J d 3 r n(r)v s [n] (r) = E s - V s [n] (4) 

where i?«(r) is the effective KS potential and E s the 
ground-state energy of the KS system, i.e., the sum of 
the KS eigenvalues. The F and T s functionals defined 
in this way can be used to establish a density variational 
principle for E v [n] and to define the exchange-correlation 
energy E xc [n] = F[n] — T s [n] — Eh[ti], where Ejj is the 
Hartree energy^ 

Thus tamed, degeneracy actually becomes helpful in 
further strengthening the foundations of DFT: on a lat- 
tice, any density can be written as a linear combination 
of densities arising from ensembles of degenerate ground 
states of a local potential, thus solving the discretized 
v — representability problem. 15 

Breakdown of mappings: Nonuniqueness. It is 
known— at least since 1983 that mapping 1, v[*B], breaks 
down if finite basis sets are used to represent the wave 
functions. Harrimani& gives both general arguments and 
an explicit example illustrating this breakdown. This 
breakdown of mapping 1 is the only one occuring already 
in the charge-density-only formulation of DFT, and it is 
manifestly an artifact of the use of a finite basis set. 

In multi-density DFTs, such as SDFT and CDFT, 
the mapping between the set of effective potentials and 
the set of ground-state densities can break down even 
in the complete basis-set limit, because inversion of 
Schrodinger's equation does not establish a unique rela- 
tion between the set of densities and the set of conjugate 
potentials. This is the so-called nonuniqueness problem 
of SDFT (and CDFT and others). Following an early 
observation of the problem by von Barth and HedinjiS 
the problem has been shown to be fundamental and per- 
vasive in recent work by Eschrig and Pickettii and by 



two of us?i& who provided explicit examples of differ- 
ent SDFT potentials sharing the same ground-state wave 
function. Ref. proposed a classification of nonunique- 
ness into systematic (arising from the existence of cer- 
tain constants of motion) and accidental (arising from 
special features of the ground state). In both cases, the 
nonuniqueness is associated with the external potential. 
Since the mapping \&[n] remains intact, and internal- 
energy functionals can be defined exclusively in terms 
of wave functions, 

F[n] = (*[n]\f+ft\9[n]), (5) 

the functionals E v [n] = F[n] + V[n], T s [n] = 
($[n]\f |$[n]>, and E xc [n] = F[n] - E H [n] - T.[n] still 
exist. 

Additional examples of both systematic and accidental 
nonuniqueness were found in CDFT— and in SDFT on 
a lattice— The extent to which nonuniqueness of the 
potentials affects various types of applications of multi- 
density DFTSj_as well as possible remedies, are discussed 
in Refs. lTl2(12l . 

Breakdown of mappings: Nonuniqueness and degener- 
acy. We have just seen that in the presence of nonunique- 
ness the mapping vfif] breaks down, whereas in the pres- 
ence of degeneracy the mapping *5f[n] breaks down. In- 
terestingly, a crucial fact has been overlooked in the 
standard analysis of either degeneracy or nonunique- 
ness: These complications can occur simultaneously! If 
a system with a degenerate ground state is treated with 
SDFT or any other formulation of DFT suffering from a 
nonuniqueness problem, none of the mappings holds, and 
no conventional HK theorem exists. In fact, it is ^f[n] 
that was used above to define F[n], T s [n] and ^^[n] in 
the absence of v[*&] (nonuniqueness), while v[^f] guaran- 
teed the existence of v[n] and thus of F[n], T s [n] and 
E xc [n] in the absence of \& [n] (degeneracy) . If both W [n] 
and i; ['I'] break down, it seems nothing is left. The break- 
down of both mappings is illustrated in Fig. ^ Three 
simple examples are given below. 

Our first example is an extension of the case of the 
noninteracting Li atomic to collinear SDFT. If the spin 
degree of freedom is included, each of the four degenerate 
states^ is additionally twofold degenerate with respect 
to S z . The set of external potentials B — 0,v = 3/r 
thus has an 8-fold degenerate ground-state manifold. 
The Slater determinants formed from the configurations 
\s 2 2p + | and ls 2 2p~ j have the same charge and spin 
densities. Again, we see that in the presence of a degener- 
ate ground state the densities do not uniquely determine 
the wave function. Differently from above, in SDFT we 
can now also consider the alternative set of external po- 
tentials B' — const ^ 0,v = 3/r. The spin-only mag- 
netic field B simply splits the ground-state manifold into 
two, one comprising the four spin-up configurations, the 
other the four spin-down configurations. The new ground 
state will be in the spin-up manifold, where the configu- 
rations ls 2 2p + | and ls 2 2p~ j remain and still yield the 
same densities. From the point of view of the mapping 
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between densities and potentials, this is simply the by 
now well knowniZii2*i2i22i2i nonuniqueness of the poten- 
tials of SDFT with respect to a weak collinear magnetic 
field. The full situation, however, is now one in which the 
densities do not determine the wave functions but only a 
(ground-state) manifold of them, and some members of 
these manifolds are ground states in more than one set 
of external potentials. The functionals \&[n], V[^] and 
V[n] thus do not exist. 

Consider next an interacting atom in an S = 1, L = 1 
state. Concrete examples are 6 C and 14 Si (with term 
3 P ) and 8 and 16 S (with term 3 P 2 ). In the set of ex- 
ternal potentials B = 0,v = Z/r the ground state of such 
systems is (2L + 1)(2S+ 1) = 9-fold degenerate. Let's de- 
note the members of this manifold as ^l z ,s z - Several of 
these, such as ^1.1 and have the same charge and 

spin densities. Hence, we have another situation in which 
these densities do not determine the wave functions but 
only the manifold. Now consider the same system in ex- 
ternal potentials B' = const ^ and v = Z/r. The 
states ^1,1, ^o,i an( i ^-1,1 remain degenerate ground 
states in this new set of potentials, and the density and 
spin density of the first and the last are still the same as 
for B' = 0. Hence, as in other examples of nonunique- 
ness, knowledge of this state alone does not determine the 
external potentials. Upon combining both observations 
we find that to a given set of ground-state densities (n, m) 
there may correspond more than one degenerate wave 
function (all in external potentials B = 0, v = Z/r), and 
all of these wave functions are also degenerate ground- 
states of the different set of external potentials (B',v). 
Again, the functionals ^[n], V[$>] and V[n] do not exist. 

Lastly, we discuss a modification of the one-electron 
example by von Barth and Hedin^fi Consider a sin- 
gle electron in the presence of an external 4-potential 
u> a p(r) = V{r)8 a p — [B(r) • cr] a/ g, where a is the vec- 
tor of Pauli matrices. Let w a p(j) be uniform along one 
spatial direction (say, x), with periodic boundary con- 
ditions along that direction separated by a distance L 
(which is topologically equivalent to confining the elec- 
tron on a ring). The two-fold degenerate ground state 
of the Hamiltonian H a p = — ^<$a/3 + Wapi?) is given 
by ^a( r ) — e ±lkx ip a (y, z), k = 2-k/L, with both ground 
states producing the same density. Furthermore, both 
ground states ^^(r) are invariant under perturbations 
w ' a 0( Y ) = V( v ){Sai3 - [m(r) ■ a/n(r)] a p}, where n(r) 
and m(r) are the ground-state density and magnetiza- 
tion, and V(r) is an arbitrary (but not too large) scalar 
potential function. Thus, we have found another case 
where both mappings (0 break down. 

Restoration of energy functionals. The question then 
arises whether the energy functionals -E [n], F[n], T s [n] 
and -Ea; C [n] can still be defined, even in the absence of all 
mappings that are conventionally considered the content 
of the HK theorem. To answer affirmatively, we consider 
two distinct cases, represented in Fig. In case I, & a 
and which both produce density n, are degenerate 
ground states in potentials V\ and V 2 . In case II only 



^b is common ground state of both potentials, whereas 
^> a is ground state only of V\ , but either an excited state 
or not even an eigenstate at all in V%. We first prove 
that case II cannot occur. In potential V\, we define the 
internal energy (not yet a functional of any density) as 

F x = E x - J d 3 rwi(r)n(r). (6) 

In potential Vi, ^> a and Vf^ are degenerate, so that E\ — 
(^ a \f + U+Vi\^ a ) = (9b\f + U + Vi\9t b ). Since both 
also produce the same density, the expectation value of 
V\ with and ^b is the same, so that 

(*„|f + L>|* a ) = (* b \f + U\^ b ) = F 1 (7) 

is independent of the choice of wave function. Next, ap- 
plying the variational principle to the Hamiltonian of sys- 
tem 2, we have 

(*a|T + U + V2\* a ) > E 2 , (8) 

with a strict inequality, as by assumption \& a in potential 
2 is not degenerate with the ground state of that system. 
From Eq. 10, and again making use of the fact that \P a 
and produce the same density, we can write this as 

(^b\f+U + V 2 \^b)=E 2 >E 2 . (9) 

The contradiction proves that case II cannot occur. This 
result is completely general, and implies the following 
theorem: Consider two degenerate ground-state wave 
functions in potential V\, f n and ^b- The constraint that 
these two wave functions have the same density guaran- 
tees that in any other potential V 2 either both are ground 
states ( and thus also degenerate ) or none of them is. We 
call this the joint- degeneracy theorem. Note that all of 
our explicit examples above respect the joint-degeneracy 
theorem. 

Even though one could formally define an F[n] func- 
tional in case II, we have just shown that that case cannot 
occur, so we only need to establish the existence of F[n] 
in case I. In analogy to Eq. 0, we define in potential V 2 

F 2 = E 2 - J d 3 rv 2 (r)n(r). (10) 

Since \l/ a and ^>b are degenerate also in potential V 2 , we 
have E 2 = (¥ a \f + U + V 2 \^ a ) = (¥&|T +U + V 3 \9 b ), 
resulting in 

F 2 = (* a \f +U\* a ) = (* b \f + tf\* b ), (11) 

where we again used that and ^>b yield the same den- 
sity. This is the same equation obtained above for F\. 
Again, this result is completely general, implying the fol- 
lowing theorem: Regardless of any possible degeneracy or 
nonuniqueness, two systems with same ground-state den- 
sity have the same internal energy F . Hence, the func- 
tional F[n] exists and is universal, i.e., independent of 
the potentials. This internal- energy theorem is consistent 
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with the constrained-search formulation of DFT, 6 which 
defines F[n] := min*_ n (*|T + ^|*), although in the 
presence of degeneracy this definition cannot be used to 
define if?[n]. 

Since the noninteracting kinetic energy T s is the inter- 
nal energy of the Kohn-Sham system, it is also a well- 
defined density functional, and E xc [n] = F[n] — Eu[n] — 
T s [n] can be constructed as usual. Finally, for a given 
external potential, the functional E v [n] then obviously 
also exists. 

Conclusions. We have shown both by general argu- 
ments and by specific examples that in the case of degen- 
eracy in multi-density DFTs all three mappings, 
V[ip] and V[n], and thus the entire body of informa- 
tion usually considered the content of the HK theorem, 



break down. The weaker joint- degeneracy and internal- 
energy theorems, however, still allow the definition of the 
internal-energy functional F[n], and thus also the func- 
tional T s [n), E v [n] and E xc [n}. Practical DFT, which 
assumes the existence of these functionals, is thus largely 
unaffected by the breakdown of the various mappings. 
However, we stress that we have only proven existence 
of the functionals, not their differentiability. In fact, in 
the presence of nonuniqueness, all these functionals are 
expected to display derivative discontinuities. 
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FIG. 1: (color online) Schematic illustration of the break- 
down of mappings that occurs in the presence of degeneracy 
(black) and the additional complication posed by nonunique- 
ness (red). V and n stand generically for sets of conjugate po- 
tentials and densities (e.g., Vf,vi and n|,nj. in SDFT). Large 
ovals are sets of functions, medium-size ovals collect degener- 
ate wave functions, and small ovals enclose degenerate wave 
functions that give rise to the same density. 




FIG. 2: (color online) The joint-degeneracy theorem (case I) 
and a situation excluded by it (case II) . Colored ovals enclose 
degenerate wave functions coming from a local potential. 



